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O-t- ABSTRACT 
O 

, Context. The 30 Doradus region in the Large Magellanic Cloud is one of the most 

, outstanding star forming regions of the Local Group and a primary target to study 

star formation in an environment of low metallicity. 

Aims. In order to obtain a more complete picture of the not yet consumed or dispersed 
(*~^ ' cool gas, we searched for line emission from molecular clouds that could be associated 

' with molecular hydrogen emission detected in the region. 

• Methods. We obtained a high sensitivity ^^CO J=2-l map with the 15-m SEST tele- 

. scope, complemented by pointed observations of ^^CO J=2-l and CS J=2-l. 

, Results. We report the discovery of a dense molecular cloud towards an embedded 

I young massive star at ~ 20" (~ 5 pc, at the distance of 50 kpc) northwest of R136, the 

compact massive central stellar cluster powering 30 Doradus in the LMC, that could be 
triggering star formation in the surrounding molecular clouds. We derived a molecular 
}_( ' mass of < 10* M0, a linear radius of 3 pc, as an upper limit, and a mean density of 



> 10 cm" for the cloud. The detection of CS J=2-l emission line indicates larger 
densities, ~ 10*^ cm~'^. The dense molecular cloud is associated with molecular 2.12 
fim H2 emission. We suggest that the observed molecular gas could be the remains of 
dense molecular material surviving the action of strong UV fields and winds in which 
the young massive star has formed. 

Key words, galaxies: Magellanic Clouds — ISM: HII Region — ISM: individual (30 
Doradus) — ISM: molecules 



1. Introduction 

30 Doradus in the Large Magellanic Cloud (LMC) is the brightest giant star forming region 
among all galaxies in the Local Group. It is 10 times more luminous than NGC3603, 
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Right ascension 

Fig. 1. Left: NIR 2.14 fim continuum image of 30 Doradus showing the position of IRSW-105, 
IRSW-107 (crosses) and the center of R136 (the X). The box shows the extent of the map displayed 
in the right panel. Right: ^^CO J=2-l emission integrated between 232 and 254 km s""*^ towards the 
mentioned IR sources. The black contour levels are 5, 6, 6.5, 7 and 8 K km s^^. The yellow contours 
correspond to the H2 near infrared emission. The black circle is the ^^CO J=2-l 24" beam. 

considered to be the most spectacular HII region in our Galaxy powered by a central 
compact massive cluster containing over 2000 Mq. 

30 Doradus is powered by a dense central compact cluster, R136. This cluster was 
initially thought to be a supermassive star (1000 Mq) to explain the measured UV photons 
required to ionize the HII region. Several very massive stars in R136 and other hot stars 
around (Walborn & Blades 1997) emit a very large flux of Lyman continuum photons which 
ionize most of the surrounding gas. forming a super-giant HII region. The stellar winds 
and possibly supernova explosions have eroded an asymmetric cavity, bounded on one side 
(the NW half) by denser material, seen in particular as molecular clouds (Johansson et al. 
1998). Today wc know that the cluster R136 consists of a compact dense core (0.1 pc in 
diameter), where more than 65 03-typc stars, the most luminous stars known, and more 
than 10 WR stars have been disentangled with i/ST observations (Massey & Hunter 1998). 
R136 has a stellar mass of 10^ Mq and a total luminosity of 10® Lq. The R136 cluster, 
the prototypical Super Star Cluster (SSC), is visible by naked eye at a distance of 50 kpc. 

In this sense, 30 Doradus is exceptional and the study of star formation and the molec- 
ular ambient in this HII region is crucial to understand the star forming processes around 
massive compact SSCs interacting with low metallicity molecular clouds, which are exposed 
to high UV radiation fields and strong winds from WR and luminous O stars. The physical 
conditions in this ISM resemble those that existed in the early universe and thus can shed 
important light on the primeval process of star formation. 

30 Doradus has been extensively studied in optical and near infrared (NIR) wavelengths 
(Poglitsch et al. 1995; Rubio et al. 1998). Several bright embedded infrared sources (Rubio 
et al. 1998), an H2O maser and numerous nebular micro-structures as revealed by Hubble 
Space Telescope (HST/WFPC2) images evidence an active as well as recent and ongoing 
star formation (Walborn et al. 1999 and references therein). 
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In particular, Rubio, Roth, & Garcia 1992 (RRG) reported for the first time an isolated 
IR source (source W9) located ~ 20" northwest of R136. Based on the source position in the 
J-H vs H-K color-color diagram and its spectral energy distribution, the authors proposed 
that W9 was most probably a young massive star. Ground based infrared observations with 
high sensitivity and improved spatial resolution allowed for the discovery of several new 
bright embedded IR sources associated with molecular clouds in this region (Rubio et al. 
1998). All of the RRG sources were identified and some of them were resolved into several 
components. RRG-W9 was resolved into two sources, IRSW-127, a very bright object with 
the largest IR excesses among those in 30 Doradus {Ks ~ 13.48, Ks-H ~ 2.85), and IRSW- 
105, displaced by 6". These two IR sources show spatial coincidence (within ^ 10") with 
the 0-type stars P621 and P600 classified as 03-06 V by Walborn & Blades (1997), who 
remarked that they might be very young objects on or near the zero-age main sequence. As 
Rubio et al. (1998) suggest, the IR sources and the mentioned 0-type stars likely represent 
the same star-formation envents. 

Star formation occurs in molecular clouds. These clouds are expected to be destroyed 
by the violent interaction of the newly born massive stars. The two IR sources were found 
at only 0.5 pc from the energetic massive cluster R136 and their proximity motivated us 
to study the molecular gas near the massive and luminous stellar cluster R136. 

In this work, we present ^^CO J=2-l, ^^CO J=2-l, CS J=2-l and ^^CO J=l-0 data 
towards IRSW-105 and IRSW-127. 

2. Observations 

High sensitivity observations in ^^CO J=2-l (rest frequency: 230.538 GHz) were made 
between January and July, 1998 using the 15m Swedish-ESO Submillimetre Telescope 
(SEST). The observations were done in position switch mode with a fixed reference off 
position free of CO emission. A narrow band AOS high-resolution (HRS) spectrometer 
with 2000 channels, 80 MHz bandwidth and 41.7 kHz channel separation (corresponding 
to 0.054 km s"^ for the ^^CO J=2-l hue) was used as back end. At 230 GHz the beamwidth 
and main bean efficiency of the telescope were 24" and 0.60, respectively. Calibration was 
done using the standard chopper wheel technique. The pointing accuracy, checked during 
the observations on RDor, was better than 2". Typical system temperatures were ^ 200 
K. The observations have an rms noise of 0.07 K, achieved after 4 minutes of integration in 
each position, which notably improve the sensitivity reached by Johansson et al. (1998) and 
Ott et al. (2008) in their ^^CO J=l-0 observations. Mapping was done with 10" spacings 
to produce fully sampled maps. The data were reduced using CLASS^ and linear and in 
few cases, third order polynomia, were used for baseline fitting. The spectra were smoothed 
to a velocity resolution of 0.25 km s~^. 

The observations covered the region between the two strongest CO molecular clouds 
mapped in 30 Doradus (Cloud 10 and Cloud 6, Johansson et al. 1998) where no CO 

^ CLASS is a GILDAS software for reduction and analysis of (sub)-millimeter spectroscopic 
data. GILDAS, the "Grenoble Image and Line Data Analysis Software", is a collection of software 
developed by the Observatoire de Grenoble and IRAM. 
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Fig. 2. i^CO J=2-l emission integrated over velocity intervals of ~ 2 km s ^. The contour levels 
are 0.2, 0.3, 0.4, 0.5, 0.6, 0.8 and 1 K km s'V The IR sources IRSW-105 and IRSW-127 are 
indicated by crosses. 

emission was detected by Johansson et al. (1998). However, Rubio et al. (1998) observed 
concentrations of vibrational excited H2 at 2.12 fiia in this region, clearly calling for more 
sensitive CO observations. In this study, we restrict ourselves to the observations covering 
the area near the central cluster R136. This region is shown in Figure 1 (left) with a 2.14 
/im continuum image of the 30 Doradus Nebula. Two crosses show the position of the IR 
sources and an X indicates the center position of R136. The box indicates the region that 
is analysed in the ^^CO J=2-l line. 

We also obtained ^'^CO and CS J=2-l spectra towards the position of the strongest 
i^CO J=2-l emission in the area, RA = 5'^38"^39.5", dec = -69°05'37.5" (J2000). The 
telescope beam at the CS J=2-l transition is 45" with a beam efficiency of 0.70. These 
spectra were reduced in a similar way as the ^^CO J=2-l observations. A ^^CO J=l-0 
spectrum towards this position was obtained from .lohansson et al. (1998). 

3. Results and Discussion 

Figure 1 (right) displays the velocity integrated ^^CO J=2-l emission between 232 and 
254 km s~^ covering an area of approximately 50" x 50" around IRSW-127 (W9) and 
IRSW-105, which arc indicated by crosses. The ^^CO J=2-l map likely shows two CO 
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clumps. The strongest one has its maximum very near IRSW-127 and on its SW border 
lies IRSW-105. The second CO clump is seen towards the northwest and does not show 
any IR sources. 

Previous studies of the ISM in 30 Doradus have shown that H2 structures appear inti- 
mately associated with bright IR sources, embedded O stars, and dense molecular clumps 
exposed to strong stellar radiation field (Poglitsch et al. 1995; Rubio et al. 1998). In Figure 
1 (right), we have included in yellow contours the 2.12 ^m molecular hydrogen emission 
extracted from Rubio et al. (1998). The II2 emission and the CO contours show a similar 
spatial distribution, with the IR source IRSW-127 associated with the brightest H2 and 
CO emission. 

Table 1 presents the molecular parameters of the strongest CO clump displayed in 
Figure 1 (right), which can also be appreciated in the ^ 244 km s~^ velocity channel map 
in Figure 2. Tmb is the peak brightness temperature, \ lsr the central velocity. Aw the 
line width (FWHM) and / the integrated line intensity. Errors are formal la value for the 
model of the Gaussian line shape. 

1/2 

To determine the size of the molecular clump we applied R — (^) As, where As is 
the linear grid spacing and N is the number of positions at which the cloud is detected above 
~5 CTrms- We obtain a radius of ~ 3 pc. Since the angular size of the considered structure 
is comparable to the angular resolution of the present observations, a linear radius of 3 pc 
can be derived as an upper limit for the size of this clump. We derived the CO luminosity 
Leo = 120 K km s"^ pc^, using Leo = -^co A^(As)2, where Ico = J Tmb dv. 

Table 1. ^^CO J— 2-1 molecular parameters of the strongest CO clump displayed in Figure 1 

(right). Tmb is in K, Vc and Ai; are in km s^^. Icq in K km s^^, R in pc and Leo in K km 

— 1 2 
s pc . 



T^mb 


V. 


Av 




R Leo 


0.65 ±0.05 


244.40 ±0.40 


4.40 ±0.50 


4.00 ±0.20 


3 120 



In Figure 2 we present channel maps of the velocity integrated CO J=2~l emission 
in steps of 2 km s~-^ towards IRSW-127 and IRSW-105. The two IR sources are located 
at the position of the peak and the SW border of ^^CO J=2-l emission in the ~ 244 
km s^^ velocity channel map. The spatial coincidence of IRSW-127 with the molecular 
clump peak and the large IR excess of this source, strongly suggests that IRSW-127 is 
embedded in this molecular clump. CO emission is also seen at lower and higher velocity 
components in a different spatial location as that of the ~ 244 km s^^ component. At 
lower velocities, the CO emission corresponds to an elongated structure extending from 
southeast to the center of the region. At higher velocities, the CO emission is associated 
with other molecular structure that contains the northwestern clump and others clumps 
farther away. This is evidence for a high degree of clumping in the region. 
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3.1. Physical properties 

Figure 3 shows the spectra of the ^^CO J=l-0 and 3—2-1, ^^GO and CS J~2-l transitions 
obtained towards the position of the strongest ^^CO 3—2-1 emission; RA = 5'^38™39.5'', 
dec = — 69°05'37.5" (J2000). The spectra were hanning smoothed and are presented as 
histograms for a better display. In aU the CO transitions three velocity components are 
seen, at 237 km s~^, 244 km s~^, and 249 km s~^, respectively, while the CS 3—2-1 
spectrum shows emission only at the two higher velocity components. Table 2 summarizes 
the line emission parameters as obtained from gaussian fit to each spectra. 




330 240 S50 360 

Velocity [km/s] 



Fig. 3. i^CO (1-0), ^^CO (2-1), "CO (2-1) and CS (2-1) spectra observed at RA = 5''38'"39.5', 
dec = — 69°05'37.5" (J2000). The intensity scale which is in Tmb has been adjusted to better show 
the emission line. 

We determined the line ratios of the strongest CO clump at the velocity component 
~ 244 km s~^. These ratios are presented in Table 3, which includes the intensity ratio 
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Table 2. Observed parameters from the spectra at tlie peak position RA = 5'^38™39.5'^, dec = 
-69°05'37.5" (J2000). 



Line emission 






Av 


/ 




(km s"-*-) 


(K) 


(km s"-') 


(K km s-l) 


l^CO J=l-0 


237.6 ±0.9 


0.11 ±0.04 


2.1 ±0.2 


0.30 ±0.12 




244.9 ±0.8 


0.50 ±0.06 


3.8 ±0.4 


1.82 ±0.25 




251.4 ±0.9 


0.25 ±0.05 


5.6 ±0.2 


1.33 ±0.25 


12 CO J=2-l 


237.8 ±0.9 


0.30 ±0.05 


5.7 ±0.5 


1.40 ±0.25 




244.4 ±0.7 


1.00 ±0.08 


4.0 ±0.5 


5.20 ±0.25 




249.7 ±0.8 


0.37 ±0.05 


5.5 ±0.7 


1.60 ±0.22 


"CO J=2-l 


236.7 ±0.3 


0.08 ±0.02 


0.4 ±0.3 


0.07 ±0.02 




244.4 ±0.3 


0.13 ±0.01 


3.1 ±0.6 


0.40 ±0.04 




249.8 ±0.2 


0.07 ±0.01 


2.1 ±0.6 


0.16 ±0.03 


CS J=2-l 


245.7 ±0.4 


0.03 ±0.01 


1.7 ±0.6 


0.07 ±0.02 




249.5 ±0.2 


0.07 ±0.01 


2.1 ±0.6 


0.12 ±0.03 



and the velocity integrated line intensity ratio at the peak position convolved to the same 
beam size. 

On the other hand, integrating the ^^CO J=l-0 data from Johansson et al. (1998) 
around ~ 244 km s^^, we fomid a weak emission at the position of our CO clump. 
Convolving the ^^CO J=2-l map to the J=l-0 angular resolution we obtained the 
^^CO J=2-l and J=l-0 integrated line ratio (i?2-i/i-o)- For the CO clump we obtained 

i?2-l/l-0 ~ 1-5. 

Table 3. Line ratios for the 244 km s"^ component at the peak position. The ^^CO J=2— 1 emission was 
convolved to the ^^CO J=l-0 and CS J=2-l beams, respectively. 



Peak Line Ratio 


-^^00(2-1) 
i^CO(l-O) 


^^CO{2-l) 
iaC0(2-l) 


-'^CO(2-l) 
CS{2-1) 


Intensity: T„b 


1 ±0.2 


8 ±1.2 


15 ±7 


Vel. Integ.: / T^j, dv 


1.4 ±0.3 


13 ±2 


38 ±30 



The value obtained in the isotopic line integrated ratio is in agreement with those for 
molecular gas in other star forming regions in the Magellanic Clouds (Israel ct al. 2003), and 
compatible with a low metallicity and strong radiation field. The multiple peak structure 
and the presence of wings in the spectra may indicate that the main contribution to the 
gas motion could be a turbulent velocity field (Falgarone et al. 1994; Garay et al. 2002). 

The detection of CS J=2-l towards this molecular clump implies the presence of dense 
gas. Generally, this 'high-density' tracer indicates densities ~ 10^ cm~^ with an upper limit 
of a few xlO^ cm^^ (Heikkila et al. 1999; Wang et al. 2009). The CS J=2-l line traces the 
dense interior of the clump, while the CO lines trace the external and less dense envelopes. 
The angular resolution of the CS J=2-l observations is a factor of two lower than the ^^CO 
J=2-l observation, which implies that the CS spectra covers the total 50" x 50" studied 
area. Thus the CS emission could be coming from several dense regions. The most intense 
CS component is seen at ~ 249 km and a second weaker component is at ~ 245 km 
s~^, very close to ^ 244 km s^^ where the CO main component is seen. These different 
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velocity components seen in the CS spectra could be an another evidence for dumpiness 
of the region. 

We determine the optical depth of the CO clump considering an isotopic abundance 
ratio ([i^CO/i^CO]) between 40 and 50. Heikkila et al. (1999) obtained [i^CO/i^CO] ratios 
of 38 and [HCO+ /ll^'^GO+] ratios of 35 ± 21 in 30Dor-10, and in the SMC N27 cloud, they 
obtained isotopic abundance ratios of 35 and 40-90, respectively. Recently, Wang et al. 
(2009) obtained the isotope ratio [^^C/^^C] ~ 49±5 in the star forming region N113 in the 
LMC. Assuming Tex = 10 K, that the surface filling factors and the excitation temperatures 
of the ^^CO and ^^CO emissions are equal, and using the ratio of the peak temperatures, 
we calculate the "'^^CO and "'^'^CO J=2-l optical depths for this cloud. We obtain t^^ ^ 6 
and r^^ 0.13, respectively. This result would not be significantly altered if T^x were as 
larger as 30 K. The ^^CO J=2-l line appears optically thick while the emission in the ^^CO 
J=2-l line is optically thin. The ratios obtained from our CO clump are similar to those 
obtained in N159 (Pineda et al. 2008). Different beam filling factors play a minor role on 
the line ratios, even in a clumpy medium (Gierens et al. 1992; Pineda et al. 2008). 

3.2. Molecular mass 

We have used three independent ways to estimate the molecular mass of the discovered 
molecular clump. First, we assume that the molecular clump is in virial equilibrium. 
Assuming that the clump has a spherical shape, a radius of 3 pc, which is an upper 
limit, we estimate Mvir ^ 1 x 10* M©, using the constant B = 190 for a density profile 
of p{r) (X 1/r (MacLaren et al. 1988). Secondly, we estimated the molecular clump mass 
from the observed ^^CO J=2-l luminosity using an H2 - ^^CO(l-O) conversion factor of 
X ^ 7 X 10^" cm^^ km~^ s, 2.5 times the Galactic value, as determined for this 
LMC region (Garay et al. 2002; Israel et al. 2003 and reference therein). We multiply this 
conversion factor by the integrated line ratio i?2-i/i-o = 1-5 obtained for the studied CO 
clump obtaining a molecular mass of Mqo ^ 2.7 x 10^ M©. Thirdly, we derived the mass 
using the ^'^CO J=2-l emission and by assuming LTE conditions. We assume optically 
thick ^^CO J=2-l emission and a beam filling factor of 1, which may not be completely 
true but allows us to make a first initial guess. We derived an optical depth, t^"^ ~ 0.12 
and using the fractional abundance of [H2/^^C0] ~ 1.5 x 10^ for this region (Heikkila et al. 
1999) we obtain Mlte < 1-5 xlO"* Mq. As a resuh, we derive a number density of H2 
molecules of Uvir > 2.0 x lO'^ cm"'^, nco ^ 0.5 x 10"^ cm"'^ and ulte ^ 2.8 x 10^ cm~^. 

Taking into account that the mass values are upper limits, in what follows we perform 
a rough comparison between them. The virial mass and the LTE mass are similar, being 
the LTE mass slightly larger, implying that the molecular clump may be approximately in 
virial equilibrium. The mass obtained from the observed CO luminosity is approximately 
one order of magnitude lower than the virial and LTE masses. Such difference may be 
indicative of the strong photo-dissociation that this molecular clump has suffered due 
to its location in extreme environmental conditions. The molecular clump may have a 
large H2 envelope where most of the mass is found with only a dense CO core which has 
survived from photo-dissociation (Lequcux et al. 1994). Another possibility to reconcile 
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the masses estimated in different ways is by assuming a conversion factor 13 times larger 
than the Galactic value, as obtained by Garay et al. (1993) for molecular clouds in the 30 
Doradus halo. Using this conversion factor we obtain a molecular mass from the observed 
CO luminosity of Mqo ^ 1-4 x 10"* Mq, in agreement with the LTE mass. 

4. Conclusions 

We discovered a molecular clump located ~ 20" northwest of the center of the compact 
cluster R136, where the IR source IRSW-127, a young massive star, is likely embedded. 
This source appears very strongly correlated with the molecular maximum. We derived a 
clump radius of 3 pc as an upper limit, a molecular mass of < 10"* Mq, and a density of 
> 10'^ cm~'^ from the ^^CO J=2-l emission. The i?2-i/i-o ^ l-S, higher than the average 
value of -R2-1/1-0 ~ 0.8 in a 9 arcmin beam, is consistent with optically thick, thermalized, 
dense gas (> 10"^ cm^"^) ready to form stars as pointed out by Sorai et al. (2001). 

The molecular clump lies in a region of H2 emission. Poglitsch et al. (1995) suggest 
that the clumpy H2 emission that they observed in 30 Doradus may be coming from photo 
dissociation regions (PDRs) and could be produced in dense ~ 10^ cm^'^ molecular clumps 
exposed to the stellar radiation field. The molecular clump studied in this work is exposed 
to an UV field of ~ 6000 times the average Galactic UV radiation field as determined using 
the Starburts99 code (Leitherer et al. 1999) for a 2 Myr cluster with a metallicity of 0.008 
at the distance of 20" of R136 (Guzman, private communication). Thus, the clump could 
have a large envelope of self-protected molecular gas. 

Our detection of CS J=:2-l emission confirms the presence of dense gas, 10® cm~'^, in 
the molecular clump. It is detected at 245 km s^^ , very close to the velocity component 
where the ^^CO J=2-l emission is stronger and spatially correlated with the two IR sources. 
Although the CS J=2-l emission is stronger at the ^ 249 km s~^ velocity component, the 
CO emission at this velocity does not show spatial coincidence with the IR sources. Thus, 
this CS emission could come from another dense clump where massive star formation has 
not yet occurred, or if it has occurred, the possible IR sources would not been detected by 
existing surveys because they could suffer an even larger extinction than that of IRSW-127. 
Higher angular resolution observations of this high density tracer would be needed to study 
the dense gas structure of the region. 

As mentioned above, the CS velocity structure shows evidence of dumpiness, which 
agrees with previous results for the 30 Doradus Nebula (Johansson et al. 1998; Heikkila 
et al. 1999) and is compatible with star forming activity. We find additional evidence of 
this activity: (1) the LTE mass (Mlte) is slightly larger than the virial mass (Mvir), which 
according to Kawamura et al. (1998) and Tachihara et al. (2000) could be interpreted as an 
indication of star formation, and (2) the isotopic line integrated ratio (^^CO/^'^CO) agrees 
with those found in different star forming regions in the Magellanic Clouds. 

We conclude that the discovered molecular clump could be the remain of molecular 
material which has survive the interaction of the strong winds and UV field of R136. This 
CO clump is an example of a dense molecular cloud with a large fraction of its mass in 
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H2 and not in CO as this molecule has suffered strong photo-dissociation from the UV 
radiation flux of the massive central cluster R136. 

We confirm the suggestions that possible energetic winds and strong UV radiation from 
the stellar cluster R136 could be triggering the star formation in the surrounding molecular 
clouds, and that IRSW-127 could be a newly formed young massive star embedded in a 
dense molecular clump. 
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